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ABSTRACT 

CTB 109 (G109. 1-1.0) is a Galactic supernova remnant (SNR) with a hemispherical shell morphology in 
X-rays and in the radio band. In this work we report the detection of 7-ray emission coincident with CTB 109, 
using 37 months of data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. 
We study the broadband characteristics of the remnant using a model that includes hydrodynamics, efficient 
cosmic ray acceleration, nonthermal emission and a self-consistent calculation of the X-ray thermal emission. 
We find that the observations can be successfully fit with two distinct parameter sets, one where the 7-ray 
emission is produced primarily by leptons accelerated at the SNR forward shock and the other where 7-rays 
produced by forward shock accelerated cosmic -ray ions dominate the high-energy emission. Consideration of 
thermal X-ray emission introduces a novel element to the broadband fitting process, and while it does not rule 
out either the leptonic or the hadronic scenarios, it constrains the parameter sets required by the model to fit 
the observations. Moreover, the model which best fits the thermal and nonthermal emission observations is an 
intermediate case, where both radiation from accelerated electrons and hadrons contribute almost equally to 
the 7-ray flux observed. 

Subject headings: acceleration of particles — cosmic rays — gamma rays: ISM — ISM: individual (CTB 109) 
— ISM: supernova remnants 



1. INTRODUCTION 

The origin of cosmic rays (CRs) with energies up to the 
knee of the CR energy spectrum is broadly attributed to the 
acceleration of particles at the shocks of supernova remnants 
(SNRs) in our Galaxy, as reviewed in |ReynoIds] ( |2008] l. Ob- 
servational evidence that suggests this is the case includes the 
detection of non-thermal X-ray emissi on from young shell- 
type SNRs, such as RX J1713 .7-3946 dKoyama et al.|[l997l 
Slane etal.|19"99l >, and Vela Jr. ( |Aschenbach|1998||Slane etaL] 
200 l| l, since these X-rays are believed to be synchrotron radi- 
ation from electrons accelerated to TeV energies at the SNR 
shock through diffusive shock acceleration (DSA). Observa- 
tions of 7-ray emission in the TeV range from some SNRs 
also support the scenario where particles are being a ccelerated 



to energies approaching 10 15 eV in these objects (Muraishi 
|et al.|2000[|Aharonian et al.|2004||Katagiri et al 
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ever, it has proven difficult to ascertain whether these high en- 
ergy photons result from leptonic emission mechanisms (such 
as inverse Compton or non-thermal bremsstrahlung emis- 
sion), or from relativi stic hadrons interacting with t he ambient 
medium (e.g. Ellis on et al.|2010t|Inoue et al.|2012| >. 

Cosmic rays interacting with regions of high-density mate- 
rial are expected to result in 7-ray emission from the decay 
of neutral pion s generated by prot on-proton interactions, as 
suggested by Claussen et al. ( 1997 1. Hence, studying SNRs 
that appear to be propagating into molecular clouds (MCs) is 
particularly important when searching for signatures of CR 
acceleration in remnants. GeV-TeV observations have proven 
that SNR/MC systems are likely an important class of 7- 
ra y sources, as supported by correlation studies carried out 
by |Hewitt et aL] d2009]). The Large Area Telescope (LAT), 
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onboard the Fermi Gamma-ray Telescope has allowed suc- 
cessful detections of several such systems in the GeV energy 
range such as SNRs W51 ( |Abdo et al.||2009c), G349.7 -0.5, 
CTB 37A, 3C 391, and G8.7-0.1 ( |Castro & Slane|20l0) . Re- 
cently, observations of the SNR W44 in the MeV-GeV en- 
ergy range with AGILE/GRID were combined with Fermi- 
LAT data and the characteristics of the resulting spectrum 
strongly suggest that the 7-ray emission originates fro m 7r°- 
decay of CR had rons interacting with ambient material (Giu- 
|liani et al.|2011) . On the other hand, observations of SNR RX 
J171 3.7-3946 with Fermi -LAT have been used t o argue that 
its 7-rays are produced fr om leptonic processes (Abd o et al.| 
|201 la| [Ellison et al.||2010] l. However, alternative scenarios 
where the SNR is expanding into complex surroundings, such 
as clumps of interstellar material, offer plausible explanations 
for the broadband characteristics of the remnant where the 7- 



rays result from pion decay emission (e.g. |Inoue et al.[ 2012). 

In addition to the information gathered through observa- 
tions of non-thermal radiation from SNRs, studies of the 
hydrodynamic evolution and the properties of the shocked 
medium and ejecta have yielded important information about 
CR production. The X-ray morphology of Tycho's SNR and 
SN 1006 indicates that the shock compression ratios in these 
object s have been modified by the acceleration o f cosmic ray 
ions ( |Warren et al.||2005j |Cassam-Chenai et al. 2008), while 
comparisons of the postshock plasma temperatures and shock 
velocities observed in SNRs IE 0102.2-7219 and RCW 86 
also point to a significant fraction of t heir explosion energy 
being placed in r elativistic particles (Hughes et al. 2000| 
|Helderetal.|2009) . 

The Galactic SNR CTB 109 (G109.1-1.0) w as originally 
discovered as an extended X-ray source by |Gregory &| 
Fahlman (1980), with the Einstein X-ray Observatory. The 



morphology of the source in the 0.1-4.5 keV energy range 
was identified as that of a semicircular shell, and corresponds 
well to that in the radio band, as observed in subsequent 
A49 cm observations by Hughes et al. (1981, 1984), with 



2 



Castro et al. 




the Westerbork Synthesis Ra dio Telescope (WS RT). Using 
XMM-Newton observations, Sasa ki et al.J ([2004) performed 
a detailed X-ray study of CTB 109, and round the emission 
from the remnant to be thermal in nature. 

The environment into which CTB 109 is expanding is a 
complex one, and is likely the source of the particular mor- 
phology of this remnant. There is a giant molecular cloud 
(GMC) located nearby, west of the SNR, which contains sev- 
eral H II regions, including the cl oud S152 (|Israel| 1 1 98Q|> . 



Tatematsu et al.| ( |1987 



17) 



From 12 CO and 13 CO observations 

conclude that CTB 109 is in contact with the GMC, and that 
the molecular material is responsible for the anisotropic SNR 
expansion. Since the SNR is a half shell, not only in the X-ray 
band, but also in radio continuum observations where the ef- 
fects of absorption by the circumstellar and interstellar media 
are negligible, the lack of emission to the west is attributed 
to th e SNR shock interac ting with the cloud of dense mate- 
rial. |Kothes, Uyaniker~& Yar (2002 ) compare the results of 
H I and CO observations towards CTB 109, with the parame- 



ters obtained for H II regions in the complex ( Tatematsu et al. 
1990), and derive a distance of Z>s nr = 3.0 ±0.5 kpc. St udies 



by |Durant & van Kerkwijk (2006 ) and |Tian et al. (2010) have 
found larger estimates of the distance to CTB 109 (7.5 kpc 
and 4 kpc, respectiv ely). Recently, this deb ate has been thor- 
oughly discussed in Kothes & Foster (2012), who determined 
the distance to be 3.2±0.2 kpc. 

In addition to the complexities of its interaction with its sur- 
roundings, CTB 109 is also of great interest because there is a 
magnetar located within it. Gregory & Fahlman ( 1980) iden- 
tified this compact source at approximately the center of cur- 
vature of the SNR shell, and it has since been determ ined to 
be an a nomalous X-ray pulsar (AXP), IE 2259+586 (Woods 
|et al.|2 004 and references therein). No radi o counterpart has 
been found t o the X-ray pulsar however (Kothes, Uyaniker 
|& Yar [2002) l . It has been suggested that the ejecta of SN 



explosions that result in magnetars as compact remnants of 
the progenitors should be more energeti c than 10 SI ergs, the 
canonical value for other SN expl osions (Duncan, & Thomp- 
|son|19 92; Buccianti ni et al.|2007[) . However, no evidence for 
explosion energies larger than 10 erg has been found in ob- 
servat ions of other SNRs w ith mag netars, s uch as Kes 73 and 
N49 ( |Vink & Kuiper|2006| ). |Sasaki et al.| ( |2004| ) determined 
the explosion energy of CTB 109 to be E = 0.7 x 10 51 a!3 erg, 
where = £>snr/(3 kpc), using the X-ray observations of the 
remnant and assuming instantaneous equilibratio n of ion and 
electron temperatures. Since [Durant & van KerkwTjk] ( |2006] > 
estimated the distance to CTB 109 to be Z)$nr ~ 7.5 kpc using 
observations of nearby red clump stars, they conclude that the 
explosion energy must have been ~7x 10 , and therefore 
that this remnant is an example of the "super-energetic explo- 
sions" expected from SN where magnetars form. The most 
recent distance estimates appear to rule out this scenario. 

Here we report on 7-ray observations of CTB 109 with 
the ferm/-LAT, and we model the broadband characteristics 
of the remnant using a simulation that includes hydrody- 
namic evolution, DSA, nonthermal emission modeling and a 
self-consistent calculation of the X-ray thermal emission ob- 
tained by tracking the non-equilibrium ionization state of the 



shocked plasma ( 


Ellison et al. 2012| Patnaude et al. 2009 Pat- 


naude etal.|2010 


and references therein). We describe several 



distinct scenarios for which the broadband characteristics of 
CTB 109 are well matched by the model using a distance of 
~ 3 kpc, where for some parameter sets the MeV-GeV emis- 
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FIG. 1. — Smoothed Fermi LAT count maps of front converted events in 
the range 2 to 200 GeV (units are 10 4 counts deg ), of the 1° X 1° region, 
centered on SNR CTB 109. The pixel binning is 0.01°, and the maps are 
smoothed with Gaussians of width 0.2°. White contours represent the X-ray 
emission (0.5-5.0 keV) from XMM-Newton data, as explained in Section 3. 
The black contour curves delimit region where TS> 25-36-49. The small 
blue circle indicates the position of IE 2259+586. 

sion observed originates in leptonic processes, others where 
7r°-decay emission is the dominant source, and a special sce- 
nario where a combination of leptonic and hadronic emission 
is required to fit the spectrum. We also briefly discuss the re- 
sults of simulations of CTB 109 assuming a larger distance of 
7.5 kpc, w hich requires the "super - energ etic explosion" pro- 
posed by Durant & van Kerkwijk (2006 ). The observations 
are described in Section 2, the model, fit parameters, and dis- 
cussion of the results are presented in Section 3, and the our 
concluding remarks are included in Section 4. 

2. OBSERVATIONS OF SNR CTB 109 

2.1. Fermi -LAT and CTB 109 

In this work, 37 months of data from the Fermi Gamma-ray 
Space Telescope Large Area Telescope (Ferm/-LAT) (from 
August, 2008 until September 2011) are analyzed. Only 
events belonging to the Pass 7 V6 Source class, which reduces 
the residual background rate as explained in detail in Acker- 
|mann et"aL[ ( |2012| l, have been selected for this study! The 
updated instrument response functions (IRFs) used are called 
"Pass7 version 6", which were develope d using in-flight data 
(|Rando et al.||2009[ |Ackermann et al.|[20T2]>. Additionally, 



only events coming from zenith angles smaller than 100° are 
selected in ord er to reduce the c ontribution from terrestrial 
albedo 7-rays (Abdo et a l.|2009b] l. All analyses include data 
from a circular region in the sky centered on the position of 
CTB 109, with radius 15°. 

The 7-ray data in the direction of CTB 109 are ana- 
lyzed using the Fermi Science Tools v9r23pln The maxi- 
mum likelihood fitting technique is employedTusing gtlike, 
to obtain m orphological and spectral information (Mattox 
|et al.|[T996| . The emission models used in gtlike include a 
Galactic diffuse component resulting from the interactions 
of cosmic rays with interstellar material and photons, and 
an isotropic one that accounts for the extragalactic diffuse 
and residual backgrounds. In this work, the mapcube file 

4 The Science Tools package and related documentation are distributed by 
the Fermi Science Support Center at http://fermi.gsfc.nasa.gov/ssc 
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gal_2yearp7v6_v0 .fits is used to describe the 7-ray 
emission from the Milky Way, and the isotropic component is 
modeled using the iso_p7v6 source . txt table. 

In order to understand the spatial characteristics of the 7- 
ray emission in the field of CTB 109, data in the energy 
range 2 to 200 GeV, and converted in the front section, are 
used. The 68% containment radius angle for normal incidence 
/ronf-selected photons in this energy band is < 0.3°. Galac- 
tic and isotropic backgrounds are modeled and test statistic 
(TS) maps are constructed using gttsmap to allow for detec- 
tion significance estimates, and to best evaluate the position 
and possible extent of the source. The test statistic is the log- 
arithmic ratio of the likelihood of a point source being at a 
given position in a grid, to the likelihood of the model with- 
out the additional source, 21og(L ps /L nu n). 

Figure [T] presents a count map of front converted events 
in the range 2 to 200 GeV, of 1° x 1° around CTB 109. 
The X-ray morphology of the remnant is presented as white 
contours (from XMM-Newton data, as discussed in Section 
3), and the black contour lines encircle the region where 
TS> 25-36-49, equivalent to a detection significance of 
5cr, 6cr, and la respectively. A 7-ray source is coincident 
with the position of the remnant with centroid (a_ooo,<$2000 = 
23 h 01 m 48 s ,+58°49'48"), with 95% confidence radius =4'.2. 
The significance of the detection, obtained from the evalua- 
tion of the peak of the TS map, is ~ 7.8o\ The residual test 
statistic map, built by modeling a point source at the best- 
fit centroid of emission, shows no evidence that the source 
is spatially extended, since the residual TS values do not ex- 
ceed 3.7cr within 1° of the centroid of emission. The possible 
extension of this source was further explored by comparing 
the overall likelihood obtained with gtlike for spatial models 
of the emission from CTB 109 as disk templates spanning 
several sizes (>di s k = 6'- 12'- 18' -24'), as well as a spa- 
tial template of the X-ray emitting material as observed by 
XMM-Newton (the X-ray observations are explained in de- 
tail in §2.2), with the likelihood of the point source model. 
The disk templates are moderately more successful at fitting 
the observations than the point source model (~ 40), with TS 
peaking at ~ 45 for the 12' disk. The disk template with ra- 
dius 18' is similarly successful (TS~ 43), butt the TS com- 
parison suggests the morphology of the 7-ray emission from 
CTB 109 is not consistent with that of the X-ray emission. 
Hence, we conclude that while we cannot rule out a point-like 
emission model, the morphology of the 7-ray source is possi- 
bly as extended as the SNR. However, the spatial distribution 
of emission is inconsistent with that of the thermal X-ray pic- 
ture. 

The study of the spectral energy distribution (SED) charac- 
teristics of the source associated with CTB 109, is performed 
using events converted in both front and back sections, and in 
the energy range 0.2-204.8 GeV. The lower energy bound is 
selected to avoid the rapidly changing effective area of the in- 
strument at low energies, and because of the large uncertainty 
below 0.2 GeV related to the Galactic diffuse model used. 
gtlike is used to model the flux at each energy bin and esti- 
mate, through the maximum likelihood technique, the best-fit 
parameters. Background sour ces from the 24-m onth Fermi 
LAT Second Source Catalog (Nol an et al.| |2012p|have been 
included in the model likelihood fits. The systematic uncer- 

5 The data for the 1873 sources in the Fermi LAT Second Source 
Catalog is made available by the Fermi Science Support Center at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/2yr_catalog/ 
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FIG. 2. — Fermi LAT spectral energy distribution of the source coincident 
with SNR CTB 109. Statistical uncertainties are shown as black error bars, 
and systematic errors are indicated by red bars. The solid line represent the 
best-fit powerlaw model to the data. 

tainties of the effective area, for the IRF used, are energy de- 
pendent: 10% at 100 MeV, decreasing to 5% at 560 MeV, 



and increasing to 20% at 10 GeV (|Abdo et al. 2009a 



Nolan 



|et al.|2 012 and references therein). As an addition to the sta 
tistical uncertainties associated with the likelihood approach, 
and the systematic errors related to the IRFs, the uncertainty 
of the underlying Galactic diffuse level is considered. This 
source of uncertainty is included in the evaluation of the sys- 
tematic uncertainties by artificially varying the normalization 
of the Galactic background by ±6% from the best-fit v alue at 
each energy bin, similarly to treatments presented in (Castro 
|& Slane|2010l . 

The spectrum of the 7-ray emission is shown in Figure |2j 
where the best-fit power law model to the data is also pre- 
sented. For energies below 800 MeV, and above 51.2 GeV 
only flux upper limits are determined from the data. The best- 
fit model yields a spectral index of T = 2.07 ± 0.12, and an in- 
tegrated photon flux above 100 MeV ofF>iooMev~ l.Tx 10~ 8 
photons/cm 2 /s. 

2.2. XMM-Newton Study 

CTB 109 was observed with XMM-Newton during period 
A01, when the entire remnant was covered with three point- 
ings, distributed South, North and East (ObsIDs 00575401, 
00575402, and 00575403, respectively). A detailed spa- 
tial and spectral an alysis of t hese observations was previ- 
ously presented in Sasaki et al. ( |2004[ ). We have reanal- 
ized the data in order to best adapt to the requirements of 
broadband modeling. Data from the European Photon Imag- 
ing Cameras (EPIC) MOS 1, and MOS 2 are studied using 
the XMM-Newton Science Analysis System (SAS) version 
1 1.0.0^] starting from the observational data files ( ODFs). Fil- 
tering and data reduction details follow those in Cast ro et al.| 
(2011a), and references therein. The effective exposure time 
after standard screening is 14.1 (14.5) ks for MOS 1 (MOS 2). 
The EPIC PN camera observations were heavily affected by 
cosmic ray noise and bad pixel columns, hence making the PN 
data difficult to match to the regions studied with the MOS 1/2 

6 The SAS package, and related documentation, is distributed by the 
XMM-Newton Science Operations Center at http://xmm.esac.esa.int/sas/ 
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FIG. 3. — Smoothed, merged XMM-Newton MOS 1/2 intensity map in the 
range 0.5-5.0 keV, of the three observations towards CTB 109. The white el- 
lipse indicates the region selected for spectrum extraction. The black ellipse 
represents the area associated with emission from the AXP IE 2259+586, 
which has been removed from the analysis, and the red (dashed) ellipse shows 
the app roximate extent of the X-ray emission enhancement known as the 
"Lobe" I Sasaki et al. 2004 and references therein). 

cameras which have better energy resolution. The PN obser- 
vations are not used in this study. 

The XMM-Newton EPIC MOS mosaic image of CTB 109, 
in the 0.5-5.0 keV band, is presented in Figure [3] To investi- 
gate the spectral characteristics of the SNR, we extracted the 
spectrum for the region shown as a white ellipse, which covers 
approximately 45% of the extent of the SNR, if we project its 
hemispheric X-ray morphology onto a full shell, with radius 
18' .5 (estimated from the distance between the IE 2259+586 
and the shock position towards the East). Emission from the 
AXP IE 2259+586 was not included, nor was emission from 
other point sources in the field. Corresponding effective area 
and spectral response files for the extended extraction region 
were created, and the ancillary response file obtained allowed 
for vignetting effects to be corrected. In order to account 
for the total XMM-Newton background, a two step process 
is performed which uses both corrected blank sky products 
provided by the XMM-Newton EPIC Background Working 
GrourFJ and data from a regi on outside the SNR shell, sim- 
ilar to the method detailed in |Arnaud et al.| ( [2002[ ). 

The resulting MOS 1/2 spectra were analyzed with the 
XSPEC software package (version 12.5.1), using an energy 
range of 0.4-5.0 keV, grouped with a minimum of 25 counts 
bin -1 , and using the \ 2 statistic (see Figure El. The emission 
above ~ 5 keV is dominated by the background. The model 
for the absorbed emission from an SNR in the Sedov-Taylor 
expansion phase SEDO V (NEIVERS 2.0) descr ibes well the 
spectral characteristics ( |Borkowski et al.|[2001| l. The best-fit 
parameters for this model are shown in Table [l] where the un- 

7 The blank sky data, obtained using observations of several sky posi- 
tions, is a combination of the Cosmic X-ray background and the non X- 
ray background, which varies in time and results from a combination of 
the interaction of cosmic rays with the detector, back ground due to so- 
lar protons, and electronic noise ^Carter & Read 2007 1. Blank sky data 
from the XMM-Newton EPIC Background Working Croup are available at 
http://xmm2.esac.esa.int/external/xmm_sw_cal/background/ 
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FIG. 4. — XMM-Newton EPIC MOS 1/2 background subtracted spectra of 
the region illustrated in Figure[3]of CTB 109. The best-fit TBABS*SEDOV 
model is shown as a histogram, and its parameters are presented in TablefT] 



TABLE 1 

Spectral Fit Parameters from MOS 1/2 data with sedov 
model 



Parameter Value 



Column Density, Nh (10 22 cm~ 2 ) 


67+ 01 


Mean Shock Temperature, kTg (keV) 


n -20+0.03 


Electron Temperature a , kT e (keV) 


n is+0.03 

U1S -0.03 


Ionization Timescale, n e t (x 10' 2 s cirT 3 ) 


1 97+O.O7 
1 - z '-0.05 


Absorbed X-ray Flux b , F_, s (erg cirr 2 s~') 


6.7 x 10-" 


Unabsorbed X-ray Flux c , F una b (erg cirr 2 s~') 


7.9 x 10-'° 


Reduced \ 2 Statistic 


3.4(499 dof) 



a kT s is the electron temperature immediately behind the shock. 
b F a bs is the average absorbed X-ray flux of the region, calculated 
in the 0.5-5.0 keV energy band, using the best-fit model to the 
EPIC MOS 1/2 data. 

c Fuaub is the average unabsorbed X-ray flux of the region, calcu- 
lated in the 0.5-5.0 keV energy band, using the best-fit model to 
the EPIC MOS 1/2 data. 

certainties quoted are the 90% confidence limits (1.6cr). In 
this work, the Tubingen-Boulder model for the absorption of 
X-rays in the interstellar med ium (ISM) TBABS is adopted 
( [Wilms, Allen & McCray|2"000| . 

2.3. Radio Observations 

As part of the Canadian Galactic Plane Survey (CGPS), 
SNR CTB 109 was studied w ith the Dominion R adio Astro- 
physical Observatory (DRAO, Kot hes et al.|2 006 ). The semi- 
circular morphology is consistent with previous observations 
in the radio band by Hughes et al. (1981,1984), as well as the 
spectral characteristics: index a = 0.5, where the radio flux 
is modeled as S„ oc v~ a , and total flux at 1.4 GHz of approx- 
imately 17 Jy. The radio spectral data used i n Section 4 is 
a combinat ion of the informat ion obtained by Hughe s~et al.| 
( [19841 ) and |Kothes etaLlp006l ). 

The CO line data used to estimate the location of molecular 
material around CTB 109 in Section 3.5 was also originally 
gathered by the CGPS, from observations of 13 CO (J=l-0) 
emission with the Five College Radio Astronomy Observa- 
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tory (FCRAO) at 115.3 GHz ( [Taylor et al.|2003> . 



3. ANALYSIS AND DISCUSSION 

3.1. Hydrodynamic Evolution and Broadband Emission 
Model 

The simulations presented here are calculated with CR- 
hydro-NEI, a code that self-consistently includes the SNR hy- 
drodynamics, a semi-analytic calculation of nonlinear DSA, 
the nonequilibrium ionization conditions behind the forward 
shock, the broadband continuum from radio to gamma-rays, 
and the X-ray line emission from the interaction region be- 
tween the contact discontinuity and the forward shock. The 
X-ray emission profile is calculated using the NEI informa- 
tion for the thermal radiation component along with the X- 
ray synchrotron emission from relativistic electro ns.The CR- 
hydro- NEI model is described in greater detail in Ellis on et al.| 
pOOT] ), Patnaude et al. (2009,2010), and references therein. 
This model also includes CR escape from the shock, and how 
this population of accelerated particles diffuses into density 
enhancements upstream. A full description of ho w this is im- 
plemented is included in Ellison & Bykov (201 1 [^] 

Through the evolution of the simulation, the NEI state of 
the shocked plasma, between the forward shock (FS) and the 
contact discontinuity, is calculated self-consistently at each 
time step, using the ionization structure, free electron nu mber 
density , and electron temperature. The Raymond & Smith 
( \911\ plasma emissivity code, updated by Brickhouse et al. 
( 1995 ) is then used to obtain the thermal X-ray emission from 
the system, for the resulting plasma characteristics (Patnaude 
et al. 2010). 

The electron spectrum from the acceleration process is used 
to calculate the synchrotron and nonthermal bremsstrahlung 
radiation from the system, as well as the inverse Compton 
(IC) emission assuming standard seed photons from the cos- 
mic microwave background. The model also considers the 
nonthermal hadronic emis sion process by calc ulating the pion 
decay emission spectrum (Kamae et al. 2006), resulting from 
collisions between protons accelerated to relativistic speeds at 
the FS and those in the surrounding material. 

3.2. Model Parameters 



TABLE 2 
CR-Hydro-NEI Input Parameters 



The distance of £>snr = 3.0 ±0.5 pc, estimated by Kothes, 
Uyaniker & Yar ( 2002) 1 from CO and H I observations, is 



adopted, and only values within the uncertainty range are con- 
sidered. Given the apparent angular size of the remnant from 
the X-ray image, the radius of the semicircle is estimated to 
be R SNR = 18' ± 1' = (16 ± l)d 3 pc, where d 3 = £>snr/(3 kpc), 

8 We note that recently an e xtensive generalization of the C R-hydro-NEI 
code has been published (e.g., Lee, Ellison & Nagataki 20121. The gener- 
alized version includes a number ot effects that are not implemented in the 
version used here. These include (1) an explicit calculation of the upstream 
precursor structure, (2) a momentum and space dependent CR diffusion coef- 
ficient, (3) an explicit calculation of MFA replacing the ad hoc parameteriza- 
tion of MFA used here, and (4) a finite Alfven speed for the particle scattering 
centers based on the amplified magnetic field. These generalizations are im- 
portant for modeling certain plasma physical processes associated with NL 
DSA, but they involve a number of new parameters that are not yet strongly 
constrained by ob servations such as those available for CTB109. |Lee, El-| 
|lison & Nagataki| pOT2) showed for SNR RX J1713, where high-quality, 
broadband data exists, that models using the generalized CR-hydro-NEI code 
give essentially idenical fits to the broadband continuum and X-ray emission 
line data as the earlier version of CR-hydro-NEI used here with a modest 
change in parameters. Except for this modest change in parameters, all of the 
conclusions reached here are unchanged in the generalized version of CR- 
hydro-NEI and we have chosen to use the simpler model described in Ellison 
et al. (2012) for modeling CTB109. 



Parameter 



Value 



SN Explosion Energy, £$N (10 51 erg) 1 

Total Ejecta Mass, (M Q ) 5 

Pre-shock Proton Temperature, Tq (10 3 K) 10 

Pre-shock Magnetic Field , B (fiG) 4.5 

Temperature of CMB Photon field , Tqmb (K) 2.725 



similar to the analysis in Sasaki et al. (2004). As previously 
mentioned, the X-ray and radio morphology of CTB 109 and 
its apparent vicinity to the GMC suggest the remnant ex- 
panded into the dense molecular material towards the west 
and rapidly decelerated in that direction, resulting in its semi- 
circular sha pe. This sce nario is consistent with the results 
obtained by Wang et al. ( 1992| l using hydrodynamic model- 



ing of this remnant evolving at the interface of the diffuse 
interstellar medium and a dense (« c ioud = 36 cm" 3 ) molecular 
cloud region. Since CR-Hydro-NEI is a spherically symmet- 
ric model, the simulations have been allowed to evolve to a 
given radius (dictated by the angular size, and the distance) 
and the resulting model fluxes for fitting have been scaled to 
appropriately account for the observed emission originating 
in only a fraction of the simulated volume. Hence, the age 
of the remnant, fsNR, is varied in the model for the different 
parameter sets to match the radius, /?snr- 

While the likely association with AXP IE 2259+586 
strongly suggests SNR CTB 109 is the result of a core- 
collapse supernova (SN), we follow the majority of work on 
this remnant and explore cases where the ambient density pro- 
file is that of a uniform circumstellar environment with con- 
stant p re-shock proton number densi ty «q, and magnetic field 
Bo- Kothes, Uyaniker & Yar (2002) found that the structure 
of the H I in the region of CTB 109 indicates that the SNR 
is not expanding into a stellar wind density profile, and hence 
suggest that the progenitor was a star of type B2/3. Addition- 
ally, the large swept up mass calculated from th e X-ray ob- 
servations of the remnant (M sw ^100 M , |Sasaki et al.|2004| 
indicates that the emission and evolution of the SNR must be 
dominated by material outside any possible stellar wind from 
the progenitor. 

A power-law ejecta density profile is assumed, as expected 
for core-collapse supernovae, with index 7 and ejecta mass, 
M e j = 5 Mq. Small variations of the characteristics of the 
ejecta, i.e., different ejecta mass values and power-law dis- 
tribution indexes, were found to have negligible impact on 
the results of the simulations presented here. This is likely 
an effect of the advanced evolutionary stage of the modeled 
SNR, when the overall properties of the SNR are governed by 
the interaction between the shock and the CSM, rather than 



the ch aracteristics of the ejecta. The analysis by Sasaki et al. 
(2004) suggests the swept-up mass of circumstellar material 
is large, as mentioned previously, and hence the contribution 
from the reverse shocked ejecta to the X-ray flux is likely 
modest and is not included in the models. 

Several other parameters characterizing the SN progenitor 
and the surrounding environment are required, including the 
SN explosion energy, £sn and the average temperature of the 
pre-shock circumstellar medium (CSM), To. Since magnetic 
field amplification (MFA) appears to be a key effect of parti- 
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cle ac celeration in young SNR shocks (e.g., Uchiyama et al. 
2007 1, an ad hoc magnification factor, B amp , is included as an 



input parameter in the model. Some of the key model param- 
eters common to all scenarios studied are presented in Table 
[2j while the age, number density of ambient protons, rela- 
tivistic electron to relativistic proton ratio, magnetic amplifi- 
cation factor, and acceleration efficiency, ecu, have been var- 
ied in order for the models to appropriately fit the observed 
radius of the remnant, the broadband nonthermal spectrum 
(both in the radio band and 7-rays), and the thermal X-ray 
emission. Detailed descriptions and discussions of the input 
parameters required for the CR-Hydro-NEI mode l are pro- 
vided in Prison et IIJ ( |2lXV7) , |Patnaude et ak1 ( [2009l ), and ref- 
erences therein. This work focuses on several sets of parame- 
ters: those where the 7-rays originate in pion decay emission 
from cosmic ray ions produced in the SNR interacting with 
the ambient medium (the "hadronic" case); those for which 
the Fermi-LPR spectrum is the result of inverse Compton 
emission from accelerated electrons (the "leptonic" scenario); 
and mixed scenarios, where the leptonic and hadronic mech- 
anisms both significantly contribute to the 7-ray flux. While 
the dominant source of 7-rays in the "leptonic" case is emis- 
sion from accelerated electrons, all models place most of the 
accelerated particle energy in protons. 

3.3. Uniform Circumstellar Medium Density Model 

For the hadronic scenario the relativistic electron to rel- 
ativistic proton ratio was set to K ep = 0.01, which is con- 
sistent with observations of cosmic ray abundances at Earth 
(Hillas 2005). No combination of ambient proton density 
and age was found to simultaneously fit the measured radius 
^snr = (16.0 ± Q.5)dj, pc at di = 1, and the broadband spec- 
trum observed. Instead, the hadronic scenario discussed here 
was obtained using «o = 1.1 cm" 3 , and ?snr = 15000 years, 
which results in an FS radius /?fs — 14.2 pc. In order to match 
the observed angular size of CTB 109 the distance must hence 
be Dsnr = 2.7 kpc in this model. In such case B- dmp = 7.0, 
and £cr = 23% are required to fit the radio and 7-ray spectral 
data. The resulting FS speed is Vps — 370 km s -1 , the total 
CSM mass swept up by the half shell is approximately 470 
Mq, the overall compression ratio is R tot ~ 4.4, the magnetic 
field immediately behind the shock is B2 — 67 /iG and the to- 
tal fraction of the SN explosion energy deposited in CR ions 
~ 0.40. The broadband nonthermal spectrum from this model 
is compared to the observations in Figure [5] (top left panel), 
and the parameters are summarized in Table 15] The total X- 
ray flux in the range 0.5-5.0 keV, estimated from the SEDOV 
model fit to the XMM-Newton MOS 1/2 data, is also shown 
in the figure. The ratio of the pion decay flux to the inverse 
Compton flux in the energy range 100 MeV to 1 TeV in this 
scenario is F^/Fic — 14. 

The observed X-ray spectrum from CTB 109 appears to be 
almost completely thermal in nature, as shown in Section 2.2, 
and the expected synchrotron flux in the X-ray band from the 
model adapts well to this fact. The thermal X-ray emission 
model obtained from the plasma emissivity code was com- 
bined with the non-thermal emission in the 0.4-5.0 keV range 
to obtain the complete X-ray emission profile. Figure B](bof- 
tom left panel) compares the resulting spectrum model,Tolded 
through the instrument responses of the MOS 1/2 instruments, 
to the XMM-Newton observations. The TBABS prescription of 
photon absorption by the ISM has been used as a multiplica- 
tive model on XSPEC as the only free component of the model 



fit to the data (A^h— 0.85 x 10 22 cirT 2 ). 

As shown in Figure 5 (bottom left panel), the computed 
thermal X-ray emission from the CR-Hydro-NEI model re- 
produces the observed X-ray spectrum, with some notable and 
important exceptions: (1) the modeled X-ray spectrum does 
not include the clearly visible emission feature at ~1.2 keV; 
and (2) the Ly-a emission from Mg, Si, and S appears to be 
overpredicted relative to the data. The model does, however, 
reproduce the soft-band emission below 900 eV as well as the 
continuum shape. Additionally, the modeled flux is consistent 
with that which was modeled in the SEDOV model. 

The feature at 1.2 keV is likely associated with Fe-L shell 
emission. At these energies, the atomic data for iron (specif- 
ically Fe XVII - Fe XXIV) i n the emissivity code are incom- 
plete ([Patnaude et al.|20"l0 1. To compensate for this feature, 
we include in the fits a Gaussian component at this energy. 
The large predicted ratio of the Lyman-a emission to the He- 
like states is endemic to the model: the hadronic model re- 
quires a high density in order to reproduce the observed GeV- 
TeV emission. In the swept-up shocked material, this results 
in more collisions per ion per unit time, which eventually 
leads to higher overall charge states. The models overpre- 
dict the flux from Lya lines relative to the triplet states in the 
hadronic model shown in Figure 5 (bottom left panel). 

Thus we conclude that hadronic processes resulting from 
the acceleration of particles in CTB 109 are a possible ori- 
gin of the 7-ray emission observed with the Fermi-LPR, but a 
distance to the SNR on the lower end of the latest estimates 
would be required. 

In the top middle panel of Figure 5, we show the best fit 
leptonic model. In order to find a leptonic dominated scenario 
which fit the broadband spectral properties as well as the the 
size of CTB 109, the relativistic electron to relativistic proton 
ratio was required to be K ep = 0.02. We no te that this value is 
consistent with that found by |Ellison et al.| l |2010) for SNR RX 
J1713. 7-3946. In such case, ?snr = 8500 years, «o = 0.2 cm -3 , 
fi amp = 2.7, and £cr = 33% are required to fit the radio and 
7-ray spectral data, and the observed radius. At the end of the 
simulation the radius of the FS is /?fs — 15.6 pc, the FS speed 
is Vfs — 780 km s , the total CSM mass swept up by the 
half shell is approximately 110 M , the overall compression 
ratio is R tot ~ 4.7, the magnetic field immediately behind the 
shock is Z?2 — 26 /iG and the total fraction of the SN explosion 
energy deposited in CR ions ~ 0.42, and F w /Fi C — 0.35. The 
resulting broadband nonthermal spectrum is compared to the 
observations in Figure [5](fo/7 middle panel). 

The expected synchrotron flux in the X-ray band from the 
leptonic dominated model is also weak in comparison to the 
XMM-Newton observation, and the thermal X-ray emission 
model obtained from the plasma emissivity code, shown in 
Figure [5] (bottom middle panel), fits the observed spectrum 
reasonably well. Much like the hadronic scenario, the ther- 
mal X-ray emission fits the data qualitatively (the absorbed 
flux in the range 0.5-5.0 keV predicted 6.7 x 10 _11 erg cm" 2 
s -1 , also matches the flux extracted with the SEDOV model), 
with notable exceptions again related to the flux ratio of He- 
like to Lya emission. In the leptonic scenario, the emission 
from Lya states is underpredicted. The reason here is that in 
the lower density leptonic scenario, the shocked plasma is not 
ionized to as high a degree as in the hadronic model. This 
model provides a reasonable overall description of the ob- 
served nonthermal broadband spectrum, thermal X-ray emis- 
sion, and radial extent of CTB 109. Thus, nonthermal ra- 
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TABLE 3 
CR-Hydro-NEI Parameters 









no 


Dsns. 


'SNR 






Rfs 




B 2 a 






Model 


Hep 


s 


(cm- 3 ) 


(kpc) 


(10 3 yr) 


^amp 


(%) 


(PC) 


Kim 


(MG) 


^cr/^sn 


F^/F lc b 


HI (Hadronic) 


0.01 





1.1 


2.7 


15 


7.0 


23 


14.2 


4.4 


67 


0.40 


14 


LI (Leptonic) 


0.02 





0.2 


3.0 


8.5 


2.7 


33 


15.6 


4.7 


26 


0.42 


0.35 


Ml (Mixed) 


0.015 





0.5 


2.8 


11 


2.6 


35 


14.6 


4.7 


26 


0.49 


1.0 



a 62 is the magnetic field immediately behind the shock. 

b F n I F\c is the ratio of the pion decay flux to the inverse Compton flux in the energy range 1 00 MeV to 1 TeV. 



diation from shock accelerated leptons may be the dominant 
source of the 7-ray emission observed in the direction of SNR 
CTB 109. 

Finally, we consider a "mixed" scenario, shown in the top 
right panel of Figure 5. The parameters are shown in Table[5J 
and the broadband and X-ray spectral models illustrated in the 
middle top and middle bottom panels respectively. The key 
parameters in this model are the K ep = 0.015, ?snr = 11,000 
years, «o = 0.5 cm" 3 , B amp = 2.6, and £cr = 35%. This in- 
put set results in a radius /?fs — 14.6 pc, FS speed Vps — 520 
km s" 1 , total CSM mass swept up by the half shell M sw « 220 
Mq, overall compression ratio /? tot ~ 4.7, magnetic field im- 
mediately behind the shock B2 — 26 /iG and Eqr/E^ — 0.49. 
In terms of the thermal X-ray emission, this model appears to 
provide the best fit. The relative flux in the modeled He-like 
and Lya lines is consistent with the observations (in contrast 
to both the hadronic and leptonic scenarios), and the modeled 
emission around ~ 1 keV also agrees well with the data, both 
in overall normalization as well as which lines are present (no- 
tably, the shape of the spectrum around 1 keV in the leptonic 
scenario does not agree with the data, while in the hadronic 
model, where the shape is consistent with the data, the model 
consistently underpredicts the flux). 

While all the models fit the broadband nonthermal spec- 
trum, the inclusion of the thermal X-ray emission allows us to 
rule out certain models. In particular, the thermal X-ray emis- 
sion favors a model where F^/Fiq. sa 1. In this case, we find 
that we can reproduce both the nonthermal emission as well as 
the shape and bulk properties o f the thermal X-ray emission. 

Since |Kofhes & Foster (2012 ) considered all other previous 
distance estimates and, using new data analysis, concluded 
that the remnant is in the Perseus Arm, at a distance of ap- 
proximately 3 kpc, we do not focus on a full treatment o f the 
£sn = 7 x 10 51 erg, £>snr = ' kP c scenario proposed by 



Du- 



rant & van Kerkwijk| ( |2006| l. However, we did explore this set 
of parameters to establish whether new constraints on the dis- 
tance could be determined by modeling the broadband char- 
acteristics. We found that we could not rule out the larger- 
distance parameter set, and that the main constraint that can be 
drawn is that the gamma-ray emission in such scenario must 
be the result of inverse Compton emission, since hadronic and 
mixed cases cannot fit the SNR characteristics appropriately. 

3.4. Other possible sources of "/-rays 

There are some alternative scenarios explaining the 7-ray 
emission detected that should be considered. One possibility 
is that the 7-rays result of high energy radiation from AXP 
IE 2259+586. Another scenario is that where the particles 
accelerated at the SNR shock are interacting with the dense 



molecular material nearby. 

The AXP IE 2259+586 has been car efully studied in the 
X-ray band since it was first detected by Gregory & Fahlman 
(1980). Observations with the Ginga satellite found a pulsar 
period of P ~ 7 s, and spin-down rate P ~ (3-6) x 10 13 s s" 1 
( jlwasawa, Koyama & Halpern|1992| l. In 2002, multiple out- 
bursts similar to those of soft-gamma repeaters (SGRs) were 
detected from IE 2259+586, and changes in the pulsed flux, 
persistent flux, energy spectrum, pulse profile, and spin-down 
of the source were also evidenced. 

The Fermi-LAT has detected a large population of 7-ray 
pulsars ( Abdo et al.||2010at |Nolan et al.|2012) , yet a recent 
search for 7-ray emission from magnetars with this instru- 
ment was unsuccessful (Abd o~et al.||2010b) l. These authors 
found no definitive evidence of a source coincident with IE 
2259+586 and they derived a flux upper limit in the 0.1-10 
GeV range of 1.7x10"" erg s" 1 cm" 2 , with TS value 15.6. 
We do not find a 7-ray source spatially coincident with IE 
2259+586 either. The local maximum of the TS is not coinci- 
dent with the magnetar position, and it is located at approxi- 
mately 0.1° SE of the AXP. Hence, while the AXP cannot be 
ruled out as the source of the 7-ray flux observed, it appears 
unlikely due to the spatial characteristics of the emission. 

Since the SNR is embedded in a complex region with large 
areas of dense molecular material, these clouds represent an 
effective target for interacting with CRs accelerated at the 
SNR shock and generating 7-rays. The two regions that ap- 
pear closest to the remnant are the large cloud which halted 
the expansion of the SNR on the west, and the cloud associ- 
ated with the X-ray bright feature called the "Lobe" (Sasaki 
et al.| [2004 and references therein), whose position is indi- 



cated in Figure [3] 

The scenario in which the 7-rays result from CRs interact- 
ing with the larger cloud to the west of CTB 109 is clearly 
ruled out by the spatial extent and location of the 7-ray emis- 
sion detected, which is concentrated on the eastern part of 
CTB 109. The relative positions of the western molecular 
cloud and the 7-ray emitting region are shown in Figure |6| 
which combines the XMM-Newton intensity map with Spitzer 
MIPS 24/im data, and CO line emission contours from ob- 
servations of the Can adian Galactic Plane Survey (Kothes, 
|Uyaniker & Yar||2002| . In this figure the infrared and CO 
line emission (at velocity -5 1 km/s) suggest that the GMC is 
mostly confined to the western part of the field, while the 7- 
ray emission is located east of the AXP (the 95% confidence 
radius of the centroid is shown as a dashed red circle). 

A fraction of the CRs accelerated at the shock of an SNR are 
expected to escape and diffuse away, and their interaction with 



nearby regions of dense material will produce 7-rays (Ellison 
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FIG. 5. — The top row shows broadband fits to radio (open squares, Kothes et al. 2006), and Fenni-LAT (black circles) observations of CTB 109 with the 
hadronic (top left), leptonic (top middle), and mixed (top right) models, all with uniform CSM density profiles. The blue triangle and horizontal bar represent the 
total X-ray flux in the 0.5-5.0 keV range estimated from the SEDOV model fit to the XMM-Newton MOS 1/2 data. The solid black line is the total nonthermal 
emission predicted by the model. The modeled spectra from inverse Compton emission (dashed), 7r°-decay (dotted), and nonthermal bremsstrahlung (dot-dashed), 
are also shown. The bottom row presents the XMM-Newton EPIC MOS 1/2 background subtracted spectra of the region illustrated in FigurefSlof CTB 109 as 
fit by the model. The histograms in the bottom left, bottom middle, and bottom right panels represent the model X-ray spectral energy distributions from the 
hadronic, leptonic, and mixed scenarios respectively. See the text for discussion and Tables[2]and[3]for model parameters. 




FIG. 6. — Combined multi-wavelength view of the region around CTB 109. 
The smoothed, merged XMM-Newton MOS 1/2 intensity map in the range 
0.5-5.0 keV is shown in cyan, while the Spitzer MIPS 24/im image is shown 
in red. The white (blue) contours correspond to the CO line emission as 
observed by the CGPS, at velocity -51 (-54) km/s, and the levels correspond 
to 0.5 and 1.5 K. The 95% confidence radius of the centroid of the Fermi-hAT 
source is shown as a dashed red circle 



|& Bykov|201 1[ and references therein). Since no significant 
emission from the western cloud is detected by the Fermi- 
LAT, we considered what constraints this places on the nature 
of particle escape from the acceleration site. For this purpose 
we used CR-Hydro-NEI and conservative estimates of the val- 
ues for the mass of the cloud (M c i ou( j ~ 10 3 M Q ), the density 
of the cloud materia l (w c ioud ~ 150 cm" 3 ), its di stance to the 
shock (/ ~ 0.2 pc) ( |Tatematsu et al.|1987| fT99Q] >, and the dif- 
fusion coefficient for CRs in this region (Dism ' 



10 26 cm 2 s - 



Fujita et al. 2010[ ). For each of the three parameter sets used in 
Section 3.3 the simulation results in ~ 3-5% of the explosion 
energy being placed in escaping particles. Based on this, we 
estimate that the resulting 7-ray flux from the Western cloud 
would hence be < 20% of the observed flux from CTB 109, 
which is below the sensitivity level of the ferm/-LAT. We con- 
clude that a significantly larger population of escaping cosmic 
rays, or an extremely high mass for the western cloud, would 
be required to produce a detectable 7-ray signature from this 

re gion. 

Sasa ki et al.| (|2006 ) used CO and X-ray observations with 
the Chandra X-ray Telescope to observe the Lobe, and con- 
cluded that there was clear evidence that the SNR shock was 
interacting with the molecular cloud complex at that position. 
These authors found three CO clouds near the Lobe with ra- 
dial velocities between -52 and -57 km/s, which can be clearly 
seen as blue contours on Figure [6] The brighter cloud, lo- 
cated to the NE of the X-ray Lobe, is also clearly coincident 
with a bright infrar ed region (as shown in Figu re [6]), IRAS 
23004+5841, which |Wouterloot & Brand] (fl987]i established 
has the IR colors of a star-forming region. [Sasaki et al. ( 2006 ) 
believe the X-ray Lobe is the result of the partial evaporation 
of the eastern cloud, and estimate the cloud properties through 
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analysis of the X-ray emission. Assuming the shocked cloud 
material fills a sphere of radius 3', and a distance of 3 kpc, 
they derive a density of «Lobe = 0.9 cm" 3 . 

If the 7-rays originate from the proton-proton interactions 
between cosmic rays accelerated at the shock of CTB 109 and 
the shocked material in the Lobe, one can derive the density 
of the shocked cloud by considering the 7-ray flux observed. 



Drury et al. ( 1994 1 derive the pion decay flux above 100 MeV 
to be F 100 Mev « 4.4 x 10- 7 (6»£ 5 



i51«c 



D k p C ) cm s , where 

£51 is the SN explosion energy in units of 10 51 ergs, the frac- 
tion of the explosion energy deposited in accelerated particles 
is represented by 9, the density of the surrounding material 
«o is in units of cm" 3 , and is the distance to the SNR 
in kpc. Using reasonable values for the parameters such as 
£51 =1,0 = 0.4, and Z\ pc = 3, as well as assuming that cosmic 
ray production in the SNR shock occurs isotropically, and that 
the fraction of the accelerated particles incident on the Lobe 
material is equal to the ratio of the angular area of the cloud to 
that of the entire SNR shell (presuming that the CRs trapped 
in the shock interact only with the portion of the cloud that 
has been overtaken by the shock), the density of the Lobe re- 
quired to produce the observed 7-ray flux (calculated in Sec- 
tion 2.1) is «Lobe ~ 120 cm -3 . Hence, there is a large discrep- 
ancy between the observational estimate of the density of the 
Lobe from th e X-ray emission («Lobe = 0.9 cm" 3 ) obtained by 
Sasak i et al.| ( |2006), and the dens ity suggested by combining 
the analysis of |Drury et al.| ( | 1 994] > and the observed 7-ray flux. 

While the Lobe density estimates obtained from the X-ray 
and 7-ray fluxes differ by approximately two orders of mag- 
nitude, it is not possible to rule out the Lobe as the origin of 
the emission above 100 MeV. Instabilities in the postshock 
flow upon interacting with dense ambient media may result in 
considerable clumping of the molecular material. The shocks 
propagating through clumps of dense material will then be ex- 
pected to quickly become radiative and their X-ray yields will 
become much less significant, resulting in density underesti- 
mates. The 7-ray flux from the interaction of cosmic rays with 
the shocked cloud is not expected to be affected by clump- 
ing, and hence high filling factors of clumped material may 
be sufficient to explain the ratio of 7-ray to X-ray flux. This 
scenario has been suggested to explain discrepancies between 
the X-ray and 7-ray characteristics of several SNRs, i nclud- 
ing RX J1713 .7-3946 ( |Inoue et al.|20l2) , and 3C391 ( |Castro| 
|& Slane|2010) . 



4. CONCLUSIONS 

X-ray and radio observations of CTB 109 reveal this SNR 
to be a semi-spherical shell, likely due to interaction with a 
giant molecular cloud in the west which halted the remnant's 
expansion in that direction. The X-ray spectrum of the SNR 



is dominated by thermal emission, and no signature of non 
thermal radiation from the shell in this energy band has been 
detected. Additionally, CTB 109 contains AXP IE 2259+586 
which is thought to be the compact remnant of the progenitor 
core-collapse supernova explosion. 

We report the detection with the Fermz'-LAT of a possi- 
bly extended 7-ray source coincident with SNR CTB 109, 
and use the broadband properties of the remnant to constrain 
the origin of this MeV-GeV emission. We use a model of 
SNR evolution and emission that includes hydrodynamics, ef- 
ficient cosmic ray acceleration, nonthermal emission mech- 
anisms and a self-consistent calculation of the X-ray ther- 
mal emission in order to consider the system's parameters re- 
quired for the SNR to be the origin of the LAT source. We 
find that the broadband radio, X-ray, and gamma-ray obser- 
vations of CTB 109 can be reasonably fit with two distinct 
parameter sets where the dominant emission mechanism be- 
hind the gamma-ray flux is either hadronic (pion-decay) or 
leptonic (inverse Compton) in nature. However, the inclusion 
of the thermal X-ray emission in the fitting results in stronger 
constraints on the origin of the 7-ray flux and indicates that 
an intermediate scenario, where both relativistic leptons and 
hadrons contribute significantly to the MeV-GeV emission, 
provides a much superior fit to the data than the extreme cases. 

Cosmic rays interacting with high-density media are ex- 
pected to result in 7-ray emission from relativistic ions collid- 
ing with ambient hadrons. Hence, we consider the interaction 
between CTB 109 and nearby molecular clouds as a possible 
origin of the Fermz'-LAT source, and while it is possible to rule 
out an association with the GMC located to the west because 
it is too far from the centroid of the 7-ray emission, we find 
that it is possible that the source is related to an interaction 
between the SNR and a region of shocked molecular material 
known as the X-ray "Lobe". Additionally, we find the associ- 
ation between the 7-ray source and the AXP IE 2259+586 to 
be unlikely due to spatial considerations. 
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